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We studied carefully the comment by S. G. Karshenboim, V. G. Ivanov and J. Chluba [ 1] to our paper "Two-photon approxi- 
mation in the theory of the electron recombination in hydrogen" [2] and found this comment erroneous. 
, 1. The errors begin with the second formula in [1J which reads 
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T A = W A °lf s , (1) 



where Fa is the total radiative width (i.e. total decay probability of state A), W^l^g is the sum of probabilities for the different 
decay channels A — Is. Strictly speaking, no such formula does exist. Eq. ([TJ has no sense until the type of the experiment is 
specified. The standard definition of the one-photon width of the level T^ 7 is 



ri 7 = £ w a-a" (2) 

A' (E A ,<E A ) 



where ^ are the one-photon transition probabilities to the lower levels. Eq. (O can be strictly derived from QED assuming 
that T^ 7 can be defined as an imaginary part of the lowest order electron self-energy correction AE^. The idea belongs to 
O \ Barbieri and Sucher [3] and the derivation one can find, for example, in [0] and in more detail in JBt]. 

In the laboratory experiments usually the one-photon transitions WP A , are measured and the level width is defined according 
^ ' to Eq. ©. The level widths defined in this way are commonly given in the textbooks. 

Of course there are exceptional cases, like the level 2s in hydrogen when the width is defined by the two-photon transition. 
The different situation arises when we want to define the transition rate from the level A to ground state via the cascade 
transition. This is directly connected with the recombination process and this is the question which we wanted to answer in our 
_C ■ paper [2]. In this case the decay rate (if there is only one decay channel) may not coincide with the level width T A J measured 
via detecting photons with definite frequency as in the experiments of the first type, described above. 

The difference between the two types of the experiments can be well understood on the simplest example: decay of the 3s 
7—1 ' level. To the ground state this level (neglecting all magnetic and electric multipole transitions) can decay via the unique cascade 
J> . 3s — s> 2p — > Is. In the laboratory experiment of the first type one can detect the photons with frequency lu — E^ s — E 2v and 
measure the transition rate W / 3 1 s 7 _ 2p which equals to W^_ 2p = 6.317 • 10 6 s -1 . Then the width I^ should be defined as 

r£ = W&_ 2p = 6.317- 10 6 *- 1 . (3) 
*^ This is the value which is presented in textbooks. 



If we consider the experiment of the second type (recombination) we have to evaluate the differential two-photon transition 
rate dW^]_ ls (w) and then obtain the total transition rate as 

E$ s —Eis 

wll_ ls = \ J dwtUM- (4) 

o 

This evaluation was performed in (01 with the result 

W3 2 s 7 _ ls = 6.285 ■ 10 6 s" 1 (5) 



The values defined by Eq. (0 and Eq. © do not coincide. Thus the equality (Q]l does not hold in this case. The transition 3s-ls 
via the cascade is a bit slowlier than the one-photon transition 3s-2p. This can happen due to the destructive interference with 
the "pure" two-photon contribution to dW 3 ]-is ^ n case °f 3rf — > Is two-photon emission the values of T 3 J and W^J differ 
more seriously: 



Al = Wl2-2 P = 6-469 .10 7 s-\ (6) 



Wl2-is = 5-814 -10 7 s-\ 



(7) 
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Now let us look what happens in case of the 3p level decay, expressed by the formula (33) in [2]. This formula was strongly 
criticized in [ 1 ] . In this case there is direct one-photon decay to the ground state, so in the experiment of the first type the photons 
with the frequencies wi = E^ p — E\ s or u>2 = E$ p — E 2s can be detected and the total width of the 3p level according to Eq. 
©is 

r£ = w#_ lfl + w#_ 2fl . (8) 

The decay rate W 7 ^"^ in case of 3p level is defined by Eq. (33) in J2]. This equation reads 

w&L = K-u + \wl;_ 2p + w*^. ( 9) 

The first term in Eq. (|9} is, of course, dominant and by the order of magnitude is close to ([8). Two other terms represent 
the contributions of two cascades 3p-2p-ls and 3p-2s-ls, i.e. 37 transitions. Both these contributions are quite small, since the 
two-photon links of the cascade transitions from 3p to Is are very slow. The transition rate via cascade is defined by the slowest 
link of the cascade, in particular by the two-photon links 3p-2p or 2s- Is. This result automatically arises in our line-profile QED 
approach and is reflected in Eq. ©. This approach appears to be useful for the analysis of the different types of the two-photon 
transitions ("two-photon approximation"). 

In yj] the authors discuss mainly the third term in the right-hand side of Eq. (O. In their Table I they suggest to replace it by 
the Wgp_2« anc ^ t° define this value as the partial width of the 37 decay. It is true that if the photon with frequency uj = E 3p — E 2s 
is detected, the decay to the ground state is necessarily 37 decay. However, in the recombination experiment we are interested in 
the transition rate 'dp — Is via 3-photon decay. This transition rate is given by the second and third terms in the right-hand side 
of Eq. ©. 

The same picture (in a more complicated form) occurs for the cascade decays of the levels 4s,4d discussed in [2]. The 
difference between T^ 7 and W^^ a was not recognized by the authors of |Q]|. No wonder that they found the difference between 
these quantities in their Tables I and II in OJJ] . 

In our paper H another question which is related to astrophysics was also discussed. This is a question of the "radiation 
escape" from the matter. We considered this problem in a very simplified formulation: every two-photon transition contributes 
to this escape. So let's look at Eq. (O from this point of view. Now we want to know how probable is the two-photon emission 
in this process. Analysis of Eq. (0 helps us to answer this question too. The third term in the right-hand side of Eq. ((9) means 
that with the probability WgpZ^s/Tsp (this i s branching ratio for the transition 3p-2s) the electron will arrive at the level 2s, 
where it has no choice but to go down to Is state, emitting two photons. This contribution to the radiation escape is, of course, 
taken into account in the astrophysical balance equations. The contribution of the second term in the right-hand side of Eq. (0, 
i.e. W^_} 2p was not vet discussed in the context of the radiation escape. However, this contribution should be quite small due to 

the small branching ratio W^ p _ 2p /T3 P . The same concerns the two-photon decay channels in Eq. (47) in their contributions 
are also negligible. Still it does not mean that such contributions are always negligible: the particular examples will be given in 
a subsequent paper. 

2. Another error concerns the method of regularization of the resonance denominators. There is a well-defined unambigous 
procedure first formulated in QED by F. Low H an d later applied to the cascade processes in (ZD. According to this approach, 
an infinite number of the lowest-order electron self-energy corrections should be inserted in the resonance approximation itno 
Feynman graph, corresponding to the two-photon decay (see Fig. 1). These insertions can be converted to the geometric 
progression. Then the graphs Fig. la), lb) produce the energy denominators, containing T^ 7 + T, 7 , where n denotes the 
resonance intermediate state. Summation of singular terms, corresponding to the insertions in Fig. lc) should be done in a 
different way. The singularities in Fig. lc) are not connected with the resonance and can be avoided only by introducing the 
Gell-Mann and Low [8] adiabatic parameter a and finally evaluating the limit a — > 0. It can be shown [9] that the summation 
of the singular graphs converts to the factor exp(—iAEf^/a) in the amplitude. For the probability of transition to the ground 
state this gives 



lim 

a-yO 



1. (10) 



Moreover, with this procedure the Lamb shift for the ground state also arrives in the energy denominators JH. 

The same approach can be unambigously applied to the 37, 47 etc transitions. However, this clearly formulated procedure does 
not satisfy the authors of $M • They do not refer neither to H not to lO] . Instead they refer to the paper lUOIl by S. G. Karshenboim 
and V. G. Ivanov, where the problem of regularization of the resonant terms is considered at the quantum-mechanical (QM), i.e. 
phenomenological level. This eventually leads to the crucial errors. For example, the authors of JlJ insist that the initial width 
should not enter the energy denominators. In the QED picture it is evident that it should since otherwise we have to omit for 
some unknown reasons the contribution of the graphs Fig. la). 
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In the paper ifioll the eigenmode decomposition of the Green function was used for the description of the sum over intermediate 
states. This decomposition containes the denominators with the phenomenological width as a parameter. The authors of ifioll 
attribute this parameter fully to the intermediate state r^ 7 , though, according to QED, in case of cascade transition it should be 
attributed to the sum r\ l7 + T^, where i is the initial state. The same is true for the every link of the cascade as it was explicitly 
proved with QED approach in H. 

This error demonstrates the danger of using the phenomenological QM theory instead of QED for the description of the 
complicated processes in atoms. It is especially regrettable since it can deteriorate in future the accurate astrophysical codes. 

3. We turn now to another conceptual error contained in [yj. This error concerns considerations about the validity of the 
theory developed in our paper [2], i.e. the theory of the cascade transitions in the resonance approximation. The authors of JU] 
draw the conclusion that these validity conditions are 

Tinitial 1 , F nna l .... 

< 1 and < 1, (11) 

-L intermediate ^ intermediate 

where the widths of the initial, final and intermediate states are involved. We should stress that the conditions Eq. (fTTT i contradict 
to the basic principles of QED and, in general, of Quantum Field Theory (QFT). The main tool of the QFT is the S-matrix which 
is constructed in such a way that one has to fix only initial and final states (stable or metastable) and then perform the uniquely 
defined evaluations that result in the transition probability (cross-section). No additional information about the intermediate 
states for any cascade (or not cascade) process is required. All this information is contained in the S-matrix. The widths of 
the intermediate states in the cascade process will arise automatically via the summation of the self-energy insertions, described 
above (see Fig. 1). No conditions of the type Eq. ( fTTT l are necessary. 

All the S-matrix elements are finite provided that the standard renormalization procedure is performed. The divergent integrals 
in the cascade transition matrix elements arise only due to the improper use of the QED perturbation theory and can be eliminated 
by the resummation of the self-energy graphs. In no way the properties of the intermediate states can influence the validaty of 
the general QED approach to the evaluation of the S-matrix elements for any transitions. Of course, in astrophysics the situation 
is much more complicated than for the free atoms, considered in [2], but the basic principles should not be violated in any case. 

The authors of fljj] also raise the question about the dependence of the initial state on creation process and connect this question 
with the validity of the cascade (resonance) approximation. In principle, the QED treatment of any process is absolutely well- 
defined only when this process starts with the stable (ground) state and ends up with this state too. Then no question of 
preparation can arise. Still it is possible also to consider the metastable states as initial ones. Note that in this context all the 
atomic states can be considered as metastable since the energy level widths are always smaller than the intervals between the 
levels. Of course, the formulation of the problem of the transition probability from a certain initial excited state is always an 
idealization. This formulation is possible provided that in the process of creation of the initial state the so-called nonresonant 
(NR) corrections can be treated as negligible. The nonresonant corrections were first introduced by F. Low ijfl and during the last 
decade were studied for different atomic process es 111 if -111 411 . The NR corrections can depend on the creation story. However, 
the general conclusion after the series of works ||IlHll4l| was that they are extremely small for all types of the excitation (i.e. 
creation) processes. Hence, we can safely assume that the creation story of the initial state in no way affects the cascade transition 
evaluations. 

Concluding, we can state that, contrary to the statement in [ 1], all the formulas and the results in our paper [2] are quite correct. 
As it was formulated in |2|] "The main goal of our paper is the formulation of the two-photon approximation which allows for the 
rigorous incorporation of all types of two-photon processes. This may be important for more accurate astrophysical investigations 
of the cosmic radiation background". Some applications of the results of paper [2], more close to astrophysics will be presented 
in a subsequent paper. 
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Fig. 1. Electron self-energy insertions into the initial, intermediate and final states for the two-photon 3s- Is transition. The double solid line 

represents the electron in the field of the nucleus, the wavy lines represent the emitted ot virtual photons. The electron self-energy loop 
represent the infinite chain of SE insertions in the resonance approximation. The insertions in Fig. la) lead to the arrival of the initial level 
width T;i s in the resonance energy denominators, the insertions in Fig. lb) correspond to the intermediate level width and the insertions of the 

type lc) produce the Lamb shift (pure real) for the ground state Is. 



